The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation          下载免费PDF全文

Seung Hyun Lee  In Wook Yeo  Kang‐Kun Lee  Won Sang Lee 《水文研究》2017,31(20):3580-3587

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.  相似文献   

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures          下载免费PDF全文

Wenjun Gong  Yunsheng Wang  Youjun Ning  Yonghong Luo  Yujie Wang 《国际地质力学数值与分析法杂志》2017,41(10):1267-1282

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Phosphorus fractions and sorption dynamics in the sediments of two Ca-SO4 water reservoirs in the central Argentine Andes     

Jos&#; Gabriel Le&#;n  Fernando Luis Pedrozo  Pedro F&#;lix Temporetti 《国际泥沙研究》2017,32(3)

  相似文献   

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion          下载免费PDF全文

Jiazhong Qian  Lei Ma  Hongbin Zhan  Qiankun Luo  Xiao Wang  Mu Wang 《水文研究》2016,30(11):1718-1726

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Discrete element model for crack propagation in brittle materials          下载免费PDF全文

Ba Danh Le  Georg Koval  Cyrille Chazallon 《国际地质力学数值与分析法杂志》2016,40(4):583-595

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method          下载免费PDF全文

A. Verde  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(6):865-886

An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds          下载免费PDF全文

K. Huang  A. Ghassemi 《国际地质力学数值与分析法杂志》2016,40(17):2293-2311

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2014年鲁甸M_S6.5地震前后地震活动性          下载免费PDF全文

侯金欣  王宝善 《地球物理学报》2017,60(4):1446-1456

2014年 8月3日在云南省鲁甸县发生M_S6.5地震，造成了巨大的人员伤亡和财产损失.关于鲁甸地震的起因，特别是鲁甸地震与溪洛渡水库之间的关系存在一些争议.分析中强地震前后地震活动性的变化可以为判断主震起源提供一定的依据.然而，由于鲁甸地震震中附近固定台站数目少、台站距离震中远及主震后尾波干扰等因素的影响，会造成地震目录的缺失，从而可能影响地震活动性的可靠分析.为了获得更加完备的地震目录，本文采用模板匹配方法对鲁甸地震前后小地震进行了搜索和识别.以台网目录中的541个事件作为模板，在主震之前33天到主震后5天的连续波形中识别出991个未在目录中列出的地震，使余震的完备震级从1.8降低到1.3.根据检测目录得到主震前33天内的b值为1.04±0.18，余震的b值在震后5天内由0.5左右逐渐上升至0.9附近并保持稳定；根据2009年到主震发生前的台网目录得到该区域的背景b值为0.93±0.04.即鲁甸地震前后区域地震活动b值均与背景b值接近，这与典型水库诱发的中强震具有的b值特征不相符.同时通过进一步分析台网目录，未发现溪洛渡水库蓄水前后研究区域内的地震活动性存在明显变化，这也与典型水库诱发地震区域的地震活动性特征不相符.综合以上结果，我们认为2014年鲁甸M_S6.5地震不具有典型水库诱发地震的特征.  相似文献   

不同岩溶储水结构分析与地球物理勘察   总被引：2，自引：0，他引：2  

甘伏平 《地质与勘探》2011,47(4):663-672

岩溶区储水结构是寻找地下水最为有利的靶区。在可溶岩与非可溶边界呈串珠状发育的岩溶泉、漏斗等,常指示地下岩溶强烈发育,接触界面的可溶岩为找水的有利部位。在背斜轴部,由于来自侧翼的挤压作用,灰岩、白云岩等硬脆性岩石易形成张性裂隙密集带,成为地下水的有利赋存部位。断裂、构造裂隙密集带使地下水互为沟通,成为良好的储水载体。在这...  相似文献   

基于相关系数和Fisher最优分割法的汛期分期研究   总被引：1，自引：0，他引：1  

王双银  杨筱筱 《水文》2011,(4):1-5

从径流形成原理出发,提出了采用降雨径流相关系数作为汛期分期指标的基本思路。以陕西省石头河水库1954~2007年的日降水和日流量资料为基础,采用主成分分析法对影响时段径流的因素进行了分析,筛选出了影响不同时段径流量的主要因素,计算了主要影响因素与时段径流量的相关系数,采用Fisher最优分割法将石头河水库汛期划分为汛前过渡期(4月1日~5月20日)、前汛期(5月21日~7月10日)、主汛期(7月11日~8月10日)、后汛期(8月11日~9月20日)和汛后过渡期(9月21日~10月31日),并结合区域气候特点和实际发生的洪水对成果进行了合理性论证。该方法资料要求低,计算简便,具有一定的推广应用价值。  相似文献